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Abstract. - The ratio of the shear viscosity (n) to entropy density (s) for the intermediate 
energy heavy-ion collisions has been calculated by using the Green-Kubo method in the framework 
of the quantum molecular dynamics model. The theoretical curve of n/s as a function of the 
incident energy for the head-on Au+Au collisions displays that a minimum region of n/s has been 
approached at higher incident energies, where the minimum rj/s value is about 7 times Kovtun-Son- 
Starinets (KSS) bound (l/47r). We argue that the onset of minimum n/s region at higher incident 
energies corresponds to the nuclear liquid gas phase transition in nuclear multifragmentation. 
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I. INTRODUCTION. - Intermediate energy heavy 
ion collisions have been extensively studied experimen- 
tally and theoretically for obtaining information about the 
properties of nuclear matter under a wide range of den- 
sity and temperature. One of the most important aspects 
of studying nucleus-nucleus collisions at these conditions 
focuses on multifragmentation and liquid gas phase tran- 
sition (LGPT) around a hundred MeV/nucleon [ThTT]. 
On the other hand, the ratio of shear viscosity to en- 
tropy density (rj/s) has been claimed to reach its local 
minimum at the phase transition temperature for a wide 
class of systems. Empirical observation of the tempera- 
ture or incident energy dependence of the shear viscosity 
to entropy density ratio for H2O, He and Ne2 exhibits a 
minimum in the vicinity of the critical point for phase 
transition [T2]. And a lower bound of rj/s > 1/4tt, ob- 
tained by Kovtun-Son-Starinets (KSS) for infinitely cou- 
pled super-symmetric Yang-Mills gauge theory based on 
the AdS/CFT duality conjecture, is speculated to be uni- 
versally valid . In ultra-relativistic heavy-ion colli- 
sions [T5H2"u] , the ratio of shear viscosity to entropy density 
was used for studying the quark-gluon plasma phase and 
a minimum value of rj/s close to the lower bounder was 
claimed. In contrast, the study on rj/s was very limited 
in intermediate energy heavy-ion collision. Therefore it 
is of very interesting to investigate shear viscosity or rj/s 
in the intermediate energy domain |21H25j . In our recent 
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paper [26], the relation between rj/s and incident energy 
or temperature is studied in the framework of one-body 
mean-filed theory, namely Boltzmann-Uchling-Uhlcnbeck 
(BUU) model, a gradual decreasing behavior of rj/s was 
observed in lower energy but tends saturated around a 
hundred McV/nuclcon. However no minimum was found, 
which may be due to the absence of the dynamical fluctua- 
tion and cluster formation in the BUU model. In contrast, 
some other signals of multifragmentation indicate a turn- 
ing point around a certain beam energy or temperature 
in experimental data as well as in some models, such as 
the quantum molecular dynamics (QMD) model. In this 
context, it will be very interesting to investigate the rj/s 
in the QMD model. 

In this work we use a microscopic transport model 
known as the isospin-dependent quantum dynamic 
(IQMD) model [57], to simulate Au+Au central colli- 
sions. The thermodynamic and transport properties are 
extracted from the nuclear matter located in the cen- 
tral sphere, with radius i?=3.5 fm. The generalized hot 
Thomas Fermi formula (GHTFF) [28h30] are employed to 
extract the thermodynamic properties, eg. temperature 
and entropy density. The shear viscosity is calculated by 
using the Green-Kubo relation. Furthermore the multi- 
plicity of intermediate mass fragments (IMFs) is also stud- 
ied as a signal of liquid gas phase transition [5Tj to verify 
the Green-Kubo's result. 

The paper is organized as follows. Section 2 provides a 
brief introduction of QMD model. In Section 3 we present 
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calculation and discussion. The conclusion and outlook is 
given in Section 4. 

II. QMD model. — The quantum molecular dynam- 
ics model approach is a many-body theory that describes 
heavy-ion collisions from intermediate to relativistic en- 
ergy |32j . The isospin-dependcnt quantum molecular dy- 
namics model [33JIMI is based on the QMD model, in- 
cluding the isospin degrees and Pauli blocking etc. Each 
nucleon in the colliding system is described as a Gaussian 
wave packet 
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Here and Pi(t) are the mean position and mean 
momentum, respectively, and the Gaussian width has the 
fixed value L=2.16 fm 2 for Au + Au system. The centers 
of these Gaussian wave packets propagate in coordinate 
(R) and momentum (P) space according to the classical 
equations of motion: 
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where (H) is the Hamiltonian of the system. 

The Wigner distribution function for a single nucleon 
density in phase space is given by 
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The mean field in IQMD model is: 

U{p) = C/ Sk y + L^Coul + C^Yuk + C4ymj 



(3) 
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where t/sky, J7coui ; C^Yuk, and C/ sym represents the Skyrme 
potential, the Coulomb potential, the Yukawa potential 
and the symmetry potential interaction, respectively |32j . 
The Skyrme potential is: 



C^Sky = a(p/p ) + (3(p/poY 



(5) 



where po = 0.16 fm -3 and p is the nuclear density. In 
the present work, the parameter set with a =-356 MeV, 
/3 = 303 MeV, and 7 = 7/6, is used, which corresponds to 
a soft equation of state. f/ Yuk is a long-range interaction 
(surface) potential, and takes the following form 
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exp(mrij)erfc(Vljm + rij/VAL)] (6) 

with V y = 0.0074GeV and m = 1.25/m" 1 . r VJ is the 
relative distance between two nucleons. The symmetry 
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Fig. 1: The time evolution of mean nuclear matter density (a) and 
kinetic energy density (b) in a central region defined as a sphere 
with radius R=3.5fm for the head-on Au+Au collisions. Different 
color line represents different beam energy which is illustrated in the 
inset. 



potential is U sym = 32 Pn pn Pp t z , where p n , p p , and po are 
the neutron, proton and nucleon densities, respectively; t z 
equals 1 or -1 for neutrons and protons, respectively. 

From Eq. © one obtains the matter density of coordi- 
nate space by the sum over all the nucleons, i.e. 
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The kinetic energy density in coordinate space could 
also be calculated from Eq. ([7]) 
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The time evolution of the mean nuclear density and kinetic 
energy density in a given central volume (with R = 3.5 fm) 
is shown in Fig. [T] Both the matter density and kinetic 
energy density are reaching their maxima around 20 fm/c. 
And the hot and dense nuclear matter survives for a longer 
time when the beam energy is lower. This can be easily 
understood as the nuclear matter experiences compressed 
and expanded more quickly at higher beam energy. 

With the help of coalescence mechanism, the fragment 
information can be given in IQMD. The intermediate mass 
fragment which is here defined with charge number greater 
than 3 and smaller than 1/3 of the system size is very 
important for nuclear multifragmentation. These frag- 
ments are larger than typical evaporated light particles 
and smaller than the residues and fission products, and 
they can be considered as nuclear fog. So the multiplic- 
ity of intermediate mass fragments (Mimfs) is related to 
the occurrence of liquid gas phase transition. Usually the 
Mimfs increases first as beam energy increases when the 
nuclear liquid phase is still dominant, and reaches a maxi- 
mum, then decreases when the nuclear gas phase becomes 
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Fig. 3: Time evolution of temperature (a) and entropy density (b) 
in the central region. Different color line represents different beam 
energy as illustrated in the inset. 



dominant [2TJ. Fig. [2] shows the Mjmfs which is ex- 
tracted from the final stage of the collision as a function 
of beam energy for head-on Au+Au collisions. One can 
see the turning energy is around 90 MeV/nucleon. If we 
study the thermodynamic evolution of rj/s as a function 
of beam energy, it is probably to exhibit a minimum at a 
certain beam energy, which means rj/s could also serve as 
a probe of nuclear liquid gas phase transition in heavy ion 
collisions. 

III. Calculations and Discussions. Thcrmody- 
namical properties of hot nuclear matter created in heavy 
ion collisions can be extracted by using the approach de- 
veloped by Faessler and collaborators [28rf30] . From a mi- 
croscopic picture of two interpenetrating pieces of nuclear 
matter, thermal quantities are deduced from the matter 
density and kinetic energy density. The extraction of the 
thermal properties is based on a generalized hot Thomas- 
Fermi formalism (GHTFF), more detailed information can 
be found in Ref. 

Time evolutions of temperature and entropy density are 
depicted in Fig. [3] Along the time scale of the collision, 
one can see that both values almost evolve isochronously, 
reach their maxima at about 20 fm/c. After the compres- 
sion stage the nuclear system begins to expand and some 
nucleons escape from the central region, and the central re- 
gion becomes cooled down. The entropy density decreases 
more quickly than temperature, this is due to direct ef- 
fect of the quick escape of the nucleons. As one can see in 
Fig.@]thc entropy per nucleon S/A in the central region al- 
most evolves iscntropically after the compressed stage for 
some times, the higher the incident energy the larger the 
S/A is. This can be understood by the transition of the 
energy of collective motion into the colliding system's ther- 
mal motion. And then S/A increases again and reaches 
a saturate value for all incident energies. Further more it 
is interesting that although thermal properties differ from 
each other during compression and expansion stage for 
different beam energies, the nuclear matter located in the 



central region has the same thermal properties at the end 
of evolution. 

Shear viscosity determines the strength of the energy 
momentum fluctuation of dissipative fluxes about the 
equilibrium state, which can be calculated by using Green- 
Kubo relation. The Green-Kubo formula for shear viscos- 
ity is defined by [37] 



d 6 r 



dt(nij(0,0)nij(r,t)), 
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where T is the temperature of the system, "0" represent 
the starting time when the system tends to equilibrium 
and t is the post-equilibration time, (7Ty-(0, 0)7Ty-(r,£)) is 
the shear component of the energy momentum tensor. In 
this work, the post-equilibration stage is defined as the nu- 
clear matter within the given central region has reached 
an equilibrium which can be judged by the stopping pa- 
rameter |38j . The stopping R p is defined as 



R 2 ^ Rt 



where R t 
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is the transverse 



and parallel momentum, resepectively. Time evolution of 
stopping is displayed in Fig. [5] From the figure, we observe 
that at the the initial stage the stopping is very small 
and its value decreases with the energy. Later there is 
a quick increase of the stopping due to the interplay of 
two body collision, mean field and Pauli blocking, and the 
collective motion energy turns into thermal energy and 
the potential energy during the compression stage. It is 
interesting to see there is a vibration of R p , this is caused 
by the transformation between the potential energy and 
kinetic energy. The stopping approaches a saturate value 
very close to 1 after oscillation, which means the nuclconic 
system in the central volume is very close to equilibrium in 
later stage of collisions. Different statring time ("0") when 
R p tends to I has been used in the Eq. © for calculating 
viscoisty. 
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Fig. 4: Time evolution of local entropy per nucleon generated in the 
central region. Different color line represents different beam energy 
as illustrated in the inset. 
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Fig. 7: The relaxation time as a function of incident energy for the 
head-on Au+Au collision for the central sphere. 



The expression for the energy momentum tensor is de- 
fined by iiij = Tij — r;5ijTi where the momentum tensor 



reads 
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Fig. 5: Time evolution of stopping in the central sphere. Different 
color line represents different beam energy as illustrated in the inset. 
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Fig. 6: (iTij (0, 0)iTij (r, t)) evolves with the post-equilibration time 
for the head-on Au+Au collision in the central sphere. Different 
symbols represent different beam energies and lines are the fits with 
Eq.[TJ] 
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with the momenta component p 1 , p 3 and the total energy 
p° of each nucleon provided by the IQMD model, /(r, p, t) 
is the phase space density of the particles in the system. In 
order to compute an integral, we assume that nuclcons are 
uniformly distributed inside the volume. Meanwhile, the 
volume with the radius i?=3.5fm is fixed, so the viscosity 
becomes 



(or 



(12) 



where r n represents relaxation time and can be extracted 
from the fit: 



(nj(0)nj(t)) oc exp ( — — ). 

7V 
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As shown in Fig. [SJ (^-(0, 0)nij(r, t)) is plotted ver- 
sus time for Au + Au collision at different incident en- 
ergies. The correlation function is damped exponentially 
with time and can be fitted by the Eq. (|13|) to extract the 
inverse slope which corresponds to the relaxation time. 
Fig. [7] shows that the relaxation time decreases as inci- 
dent energy, indicating that the system can approach to 
equilibration faster at higher incident energy, which is con- 
sistent with the results from the stopping parameter. 

Finally, we come by the shear viscosity and its ratio 
to entropy density. Fig. EJa) shows shear viscosity r\ as 
a function of incident energy. It is very interesting to 
see that the shear viscosity alone also exhibits a mini- 
mum near 120 MeV/nucleon. Since the temperatures in 
post-equilibration stage are almost the same for the colli- 
sions at different energies, the value of shear viscosity are 
mostly influenced by the tensor correlation, which reflects 
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Fig. 8: Shear viscosity r\ (a) and the ratio of shear viscosity to 
entropy density (b) as a function of incident energy for the head-on 
Au + Au collisions. 

the fluctuation of dissipative fluxes. Fig. [SJb) displays 
rj/s as a function incident energy for the head-on Au + 
Au collisions from 70 to 200 MeV/nucleon. The rj/s de- 
creases quickly with the incident energy up to a platform 
of the minimum value at around 120 MeV/nucleon and 
afterwards the curve weakly rise up slowly. The minimum 
value is around 0.6, i.e. about 7 times KSS bound (l/47r), 
which is very close to the saturated value in our previous 
BUU calculation [26]. As we expected, all rj/s values from 
the present QMD model are lager than the KSS bound. 
On the other hand, we noticed that the location of energy 
with the minimum Ntmfs and rj/s is not exactly the same, 
it might be due to our Green-Kubo calculation focuses on 
a given central nuclear region, and the IMFs represent the 
whole colliding system. 

IV. Conclusion. — In this Letter, we performed a 
first calculation on the ratio of shear viscosity to entropy 
density in the framework of quantum molecular dynam- 
ics model with the Green-Kubo method. For an example, 
the head-on Au+Au collisions are simulated for beam en- 
ergies from 70 to 200 MeV/nucleon. The generalized hot 
Thomas-Fermi formalism is used to extract the thermal 
properties of the given central nuclear system, and the 
shear viscosity is calculated with the Green-Kubo relation. 
The results display that a minimum rj/s region, about 7 
times Kovtun-Son-Starinets (KSS) bound (l/in), exhibits 
at around 120 MeV/nucleon, which may correspond to a 
liquid gas phase transition for a hot nuclear matter. The 
ratio of shear viscosity to entropy density is always larger 
than 1/47T, supporting the result of AdS/CFT. Finally, 
we would mention that the effect of equation of state and 
symmetry energy on the rj/s is also an interesting topic 
[40 43], which is beyond the scope of this Letter, will be 
studied in near future. 
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